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Abstract. 

We explore photometric properties of hypothetical iron core white dwarfs and compute their expected colors in UBVRI Johnson 
broadband system. Atmospheres of iron core WDs in this paper consist of pure iron covered by a pure hydrogen layer of 
an arbitrary column mass. LTE model atmospheres and theoretical spectra are calculated on the basis of Los Alamos TOPS 
opacities and the equation of state from the OPAL project, suitable for nonideal Fe and H gases. We have also computed UBVRI 
colors of the models and determined an area on the B — V vs. U — B and B - V vs. V - I planes, occupied by both pure Fe, and 
pure H model atmospheres of WD stars. Finally, we search for iron core white dwarf candidates in the available literature. 
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1. Introduction 

Chemical composition of white dwarf stars cores has not yet 
been convincingly determined. While most white dwarfs are 
believed to have pure carbon or mixed carbon-oxygen interiors, 
there are suggestions that cores of some WDs are built of heav- 
ier elements, and pure iron cores represent the extreme case. 

An idea that some white dwarfs have cores built of heav- 
ier elements was put forward by Provencal et al. (1998) and 
Provencal & Shipman (1999). In particular, they noted that 
the nearby white dwarf Procyon B has relatively normal white 
dwarf mass of 0.594 M , but has a very small radius (cf. also 
Girard et al. 2000). This would imply that Procyon B has a 
very heavy core, perhaps of iron. They found that radii of other 
WD stars: EG 50 and GD 140 are also located significantly 
below the mass-radius relation predicted for typical carbon in- 
terior compositions. Only recently, Provencal et al. (2002) on 
the basis of precise Hubble spectral measurements recognized 
that Procyon B is a normal white dwarf, belonging to the DQZ 
spectral class. 

Panei et al. (2000) have obtained thorough theoretical 
mass-radius relations for iron-dominated cores. They show that 
EG 50, GD 140 and 40 Eri B satisfy these relations for iron 
cores, contrary to mass-radius relations for carbon and oxygen 
cores. Their result convinced us that the existence of the iron 
core white dwarfs is possible. The above hypothesis is a suffi- 
cient reason for us to perform and publish this research paper. 
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Atmospheres of these stars are probably built from pure 
iron gas covered by hydrogen (and/or helium) layer, since strat- 
ified atmospheres are common among the known white dwarfs. 
We expect that most likely iron and hydrogen are not perfect 
gases, at least in cooler WDs. 



In the following sections we describe model atmosphere 
calculations in the local thermodynamic al equilibrium (LTE), 
which are based on a frequently used equation of state com- 
puted by the OPAL Project for nonideal hydrogen (Rogers, 
Swenson, & Iglesias 1996) and iron gas (Rogers 2002). Our 
calculations use exclusively the Los Alamos TOPS monochro- 
matic absorption and scattering opacities, consistent with the 
nonideal EOS for both gases (Magee et al. 1995). Model atmo- 
sphere and spectrum calculations are restricted to atmospheres 
in radiative equilibrium, corresponding to effective tempera- 
tures T eff > 20000 K. 



Los Alamos opacity tables were originally designed to 
study physics and pulsational instabilities in stellar envelopes 
by means of Rosseland mean opacities. Therefore the tables do 
not have spectral resolution adequate to study details of stellar 
atmospheric spectra. This is the reason that our paper presents 
and discusses mostly spectral indices U - B, B - V, V - R, 
and V - I of the broadband Johnson photometric system. We 
attempt to find a prescription for differentiating most common 
DA white dwarfs with pure hydrogen atmospheres from iron 
core white dwarfs with pure iron or stratified hydrogen-iron 
atmospheres. 
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2. Numerical methods 

2.1. Model atmosphere computations 

Model atmosphere equations are based on the following equa- 
tion of radiative transfer, suitable for planar geometry 

A* =I v -e v B v -{l-e v ) J v , (1) 

OTy 

where variables 7 V , J Y , B v and fi have their usual meaning. The 
three monochromatic variables are expressed in units erg cm -2 
sec -1 Hz -1 str -1 . Dimensionless absorption is e v = k v /(k v + 
o-y), where k v and <x v denote LTE true absorption and coherent 
scattering coefficients, respectively. Absorption coefficient k v 
has been corrected in our code for stimulated emission as in 
LTE. 

Computer code used here has been derived from the code 
atm21 adapted to coherent scattering in stellar atmosphere (ex- 
act description in Madej 1991, 1994, 1998; Madej & Rdzariska 
2000; Stage, Joss, & Madej 2002). 

The foregoing papers explain in detail the equations of hy- 
drostatic and radiative equilibrium, which are imposed on the 
model atmospheres, as well as two boundary conditions asso- 
ciated with the equation of transfer. The equation of transfer 
and the radiation field was solved with the method of variable 
Eddington factors (Mihalas 1978). 

Model atmosphere equations were expressed in discrete 
form on frequency grid v,-, i = 1,. . .,7 and standard optical 
depth grid Tj,d — 1 , . . . , D. Here, in our models, the monochro- 
matic optical depth at the fixed wavelength of 1500 A is used. 
For each model we set I = 749 and started computations from a 
very small standard optical depth point Ti = 10 -10 . We have as- 
sumed a value for D equal to D — 150 for pure hydrogen mod- 
els, while for pure iron and stratified H-Fe models D — 155 was 
set. Both the monochromatic absorption and scattering coeffi- 
ciens, and hence the e v , were determined at each frequency V; 
and depth points Td by direct interpolation of extensive TOPS 
opacity tables. The equation of state for nonideal gas, i.e. the 
values of gas density for a given temperature and gas pressure, 
were similarly determined by interpolation of the OPAL EOS 
tables either for H or pure Fe gas. 

2.2. Equation of state and opacity data 

The nonideal equation of state for partially ionized plasma, 
which was used in our research, has been extensively described 
by Rogers (1994), Rogers, Swenson & Iglesias (1996), cf. also 
references therein. These authors refer to their derivation as to 
the activity expansion of the grand canonical partition function 
of plasma, the latter being basic function in statistical thermo- 
dynamics. The method includes quantum mechanical two-body 
and many-body interactions, and electron degeneracy along 
with other effects. We feel not enough competent to present and 
discuss details of their method in the present paper. However, 
we mention after Rogers et al. ( 1 996) that their derivation of the 
EOS for various plasmas basically differs from common meth- 
ods which are based on minimization of free energy in plasma. 

This research uses EOS data which is available for a gen- 
eral user. EOS tables for pure hydrogen plasma were taken 



from the Lawrence Livermore National Laboratory www site 
(http://www-phys.llnl.gov /Research/OPAL/index.html). With 
EOS tables for pure iron gas we were kindly provided by Dr. 
F.J. Rogers (2002). 

Extensive TOPS monochromatic opacity tables have been 
prepared at the Los Alamos National Laboratory (Magee et al. 
1995). We have collected hydrogen and iron data from their 
web page at http://www.t4.lanl.gov/opacity/tops.html See that 
www site for more details on the TOPS opacity calculations. 

2.3. Synthetic colors of the Johnson UBVRI system 

Having calculated the models, we acquired their broadband 
Johnson U BVRI colors using a color synthesis procedure de- 
scribed by Bergeron et al. (1995); see also detailed review of 
photometric systems by Girardi et al. (2002). 

At first, stellar magnitudes of a theoretical spectrum are 
computed according to the standard formula 

oo 

/ FxSxdA 

m Sl = -2.5 log — + m s o , (2) 

/ S A dA 
o 

where we integrate the product of the flux energy spectrum F A 
of a model atmosphere with the transmission function S A of 
the respective broadband filter. The constant m s o can assume 
an arbitrary value. 

We used the Johnson UBVRI transmission functions de- 
fined by Bessell (1990), and restricted our research to analysis 
of colors, i.e. differences of the respective magnitudes. Then 
the definition of the respective colors, e.g. B - V color, is as 
follows: 

oo 

/ F A S B (A)dA 

B-V = -2.5 1og^ + C B -v, (3) 

/ F A S v (A)dA 
o 

cf. Bergeron et al. (1995). Constant Cb-v includes all quantities 
which are not related to flux F A . 

Note that plane-parallel model atmosphere codes compute 
theoretical spectra of radiation F A which are emitted by unit 
surface on the star, this is 1 cm 2 in case of the atm21 code. 
Therefore theoretical flux F A in Eq.[3]is never directly compa- 
rable with e.g. the flux measured on Earth or absolute flux in 
that filter, which would be measured at 10 pc distance. 

Eq. |3] can be applied to Vega, which is the fundamen- 
tal spectrophotometric standard. In our research each color of 
Vega is equal to zero by assumption. Setting F A to the flux 
of Vega allowed us to determine constants Cu-b, Cb-v, Cv-r 
and Cv-i for all theoretical colors presented in the following 
Sections. 

For the above calibration purposes we adapted absolute 
fluxes of Vega observed by Hayes (1985), for A > 4500 A. His 
paper gives values of the Vega flux in equidistant intervals of 25 
A and therefore gives poor description of the flux in the near ul- 
traviolet, where prominent higher Balmer lines merge and the 
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wavelength depence of Vega flux is very complex. Therefore 
for A < 4500 A we have taken theoretical fluxes from Castelli 
& Kurucz (1994). They have computed atlas12 model atmo- 
sphere and spectrum of Vega with a high spectral resolution, 
normalised in the same way as Hayes' (1985) observations. 

We have taken this approach to reference spectra of Vega 
from Bergeron et al. (1995), in order to make our color indices 
for hydrogen and iron atmospheres fully compatible with their 
tables of Johnson photometry colors computed for DA and DB 
white dwarfs. 

3. UBVRI colors of pure H atmospheres 

In this section we construct and discuss two different series 
of pure hydrogen model atmospheres of white dwarfs. This is 
done in order to compare and measure the discrepancies be- 
tween each other and assess the uncertainty of colors of atmo- 
spheres computed with the following two basically different 
numerical methods. By doing this, we also test the compati- 
bility of both OPAL and TOPS data with the existing model 
atmosphere codes. If effect, we draw a conclusion that our pure 
Fe model atmospheres, based on those data, are consistent with 
pure H atmospheres. 

For calculating models of hot WD atmospheres consisting 
of pure hydrogen, the following approaches were undertaken: 

(1) We assumed that plasma of hot hydrogen atmosphere is 
described by the EOS of ideal gas. Under this assumption, 
thorough models were calculated in the NLTE regime. Models 
and spectra computations were performed using tlusty195 
and synspec42 codes (Hubeny 1988; Hubeny & Lanz 1992, 
1995). Three spectral series of hydrogen (Lyman, Balmer, and 
Paschen) were considered as we computed outgoing flux spec- 
tra and color indices. Since Paschen lines are not present in the 
standard versions of both codes, we have adopted to the codes 
the opacities of thermally & pressure-broadened Paschen lines 
from Lemke (1997). 

(2) Alternatively, we assumed more factual, nonideal EOS for 
gas and derived a number of models in LTE. Tables for non- 
ideal EOS were taken from the OPAL project, while the Los 
Alamos TOPS monochromatic opacity tables were used. 

In the present work we are concerned about isolated hot 
WD atmospheres with effective temperatures T e ff in the range 
of 20000 to 100000 K and surface gravities ranging from 
logg = 7.0 to 10.0 (cgs units). White dwarfs of as high sur- 
face gravity as logg = 10.0 are not yet known, however, one 
can expect that iron core white dwarfs of a given mass have 
a particularly small radius and therefore high surface gravity. 
Consequently, we appended high gravity hydrogen model at- 
mospheres to obtain grid consistent with pure iron high gravity 
atmospheres. The radiative and hydrostatic equilibria are im- 
plied in all models. 

Grids of color indices for both NLTE (tlusty195) and 
LTE (Los Alamos) series of models were plotted on two-index 
plane, U - B vs. B - V, and are displayed on Figure 1. Colors 
of the TOPS-based LTE models are in a reasonable agreement 
with the corresponding colors of the NLTE models, especially 
for highest effective temperatures T e ff- The largest difference 
in the U - B index is equal to almost 0.1 mag (for T ef f = 20000 



Table 1. Broadband colors in stellar magnitudes for pure hy- 
drogen atmospheres. Models are computed also with nonideal 
EOS and Los Alamos (TOPS) opacities. Temperatures T e ff are 
given in thousands K and log g are in cgs units. 

Z T eff \ogg U-B B-V V-R V-I 

1 20 7 : 0895 11051 4L098 -0.206 

1 20 8 -0.965 -0.012 -0.103 -0.202 

1 20 9 -1.017 0.014 -0.109 -0.199 

1 20 10 -1.076 0.026 -0.112 -0.180 

1 30 7 -1.121 -0.193 -0.135 -0.290 

1 30 8 -1.157 -0.171 -0.139 -0.288 

1 30 9 -1.186 -0.163 -0.142 -0.288 

1 30 10 -1.211 -0.168 -0.145 -0.275 

1 50 7 -1.238 -0.284 -0.146 -0.325 

1 50 8 -1.244 -0.277 -0.147 -0.325 

1 50 9 -1.251 -0.275 -0.148 -0.323 

1 50 10 -1.256 -0.275 -0.149 -0.320 

1 70 7 -1.257 -0.313 -0.145 -0.338 

1 70 8 -1.262 -0.308 -0.150 -0.339 

1 70 9 -1.268 -0.304 -0.150 -0.336 

1 70 10 -1.271 -0.304 -0.152 -0.336 

1 100 7 -1.279 -0.329 -0.142 -0.348 

1 100 8 -1.280 -0.327 -0.145 -0.350 

1 100 9 -1.282 -0.324 -0.146 -0.349 

1 100 10 -1.282 -0.324 -0.150 -0.349 



Table 2. Broadband colors in stellar magnitudes for pure hy- 
drogen NLTE atmospheres (results of tlusty 195). Units are 
the same as in Table 1 . 



Z T eff 


l°g# 


U-B 


B-V 


V-I 


1 20 


7 


-0.812 


-0.091 


-0.209 


1 20 


8 


-0.866 


-0.035 


-0.210 


1 20 


9 


-0.927 


0.029 


-0.203 


1 50 


7 


-1.217 


-0.291 


-0.330 


1 50 


8 


-1.224 


-0.278 


-0.331 


1 50 


9 


-1.229 


-0.258 


-0.327 


1 100 


7 


-1.270 


-0.328 


-0.350 


1 100 


8 


-1.271 


-0.322 


-0.352 


1 100 


9 


-1.269 


-0.309 


-0.348 



K and logg = 8); in the B - V it is mere 0.034 mag (for 
T eff = 20000 K and logg = 7). The differences in V - I color 
are smaller still. 

We conclude that the models, based on TOPS and OPAL 
tables, are in satisfactory agreement with the NLTE tlusty 195 
models, as long as pure hydrogen models are concerned. 
Therefore we believe that pure iron or stratified hydrogen-iron 
models with nonideal EOS would be consistent with NLTE 
ideal gas models at the high T e ff limit. 

It is interesting to note that colors of pure H model atmo- 
spheres computed either with the EOS of ideal gas or with the 
OPAL nonideal EOS exhibit systhematic differences for lower 
effective temperatures T e ff of our grid. Only at the highest T e ff 
both series of models converge on the U - B vs. B - V plane. 
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Fig. 1. Comparison of theoretical col- 
ors, B - V vs. U-B, for two grids 
of pure hydrogen, hot DA atmospheres. 
Models are computed either with the 
tlusty195 and synspec42 NLTE codes, 
or with our atm21 code based on the 
nonideal EOS and Los Alamos opac- 
ities. Consistency of both families of 
broadband B — V and U-B indices is 
not perfect but is satisfactory. 



Such systhematic differences are not present in the B - V vs. 
V - I diagram and hence are not displayed here. 

Furthermore, our synthetic NLTE colors of pure hydrogen 
white dwarfs are nearly identical with both our LTE colors 
of these stars (computed with other option to the tlusty195 
code and still using ideal EOS) and the LTE colors obtained 
by Bergeron et al. (1995) in their intensive study on ages, lumi- 
nosities, and other photometric parameters of DA white dwarfs. 

3. 1 . Discussion of discrepancies 

The above-mentioned systhematic differences could be caused 
by the effects of deviations from the ideal EOS for hydrogen 
gas, since they vanish at the highest T e ff- In order to verify 
such a hypothesis, we have chosen several temperature T and 
gas pressure P g points across our sample pure H model atmo- 
sphere, T e ff = 20000 K and \ogg = 8.0, computed with the 
nonideal EOS and TOPS opacity data, see Tabled 

We compared two different densities corresponding to the 
same pair of T, P g : density Prog, obtained by direct interpola- 
tion of the OPAL tables, and p ID corresponding to an ideal H 
gas, obtained from the relevant subroutines of the atm21 code. 
One can see that the difference between both densities reach 
0.3 % at the bottom of the model, and it is much smaller in 
upper layers. Therefore we conclude that the EOS of an ideal 
gas is well fulfilled in all our pure H atmospheres, including 
also models of higher T e /f, of course. This firmly suggests that 
the impact of nonideal gas effects on the color discrepancies 



between our H models and these obtained with tlusty195 is 
negligible. 

The referee pointed out that such differences could be due 
to NLTE effects considered by the code tlusty195, or other 
reasons. Earlier studies on this subject have already predicted 
the nonideal EOS's and NLTE effects in hydrogen atmospheres 
of white dwarfs. Wesemael et al. (1980) in their extensive paper 
on WD model atmospheres and spectra, have also studied the 
influence of hydrogen line blanketing and NLTE effects on both 
U-B and B — V color indices. They found, for a pure H model 
atmosphere of T eff - 50000 K and logg = 8.0, that both &(U - 
B) and a(B - V) line-blanketing corrections are of the order 
0.001 mag (at most). 

Therefore NLTE corrections to UBV are just negligible 
in case of pure H white dwarfs, for relatively hot objects of 
T e ff > 20000 K. This conclusion is also supported by our own 
computations, which were not detailed in this paper. 

The influence of blanketing by Balmer and Lyman lines is 
much more important. Wesemael et al. (1980) have determined 
that both a(U - B) and a(B - V) corrections, LTE minus NLTE, 
are equal to -0.017 mag and +0.043 mag, respectively, for a 
model of T e jj = 50000 K and logg = 8.0. Both corrections 
should increase in their absolute value with decreasing T e jj (see 
p. 300 of that paper). The quantity of both corrections, when 
extrapolated down to the T e ff - 20000 K, seems comparable 
with both differences a(U - B) and a(B - V), cf. Fig. 1. 

Note that both series of pure H model atmospheres dis- 
played in Fig. 1 take into account the blanketing by Balmer 
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and Lyman lines. However, one can expect that perhaps TOPS 
hydrogen opacity tables do not include most sophisticated the- 
ory of H line broadening, thus causing a(U - B) and a(B - V) 
errors seen in Fig. 1. The difference strongly decreases in V - 1 
color, since in both V and / filters there are only single B aimer 
lines, and their broadening gets unimportant for the wide-band 
photometry. 

At this point we turn the reader's attention to the paper 
by Hubeny, Hummer & Lanz (1994), who discuss theory of 
NLTE model atmospheres with line blanketing also near the 
series limit (the occupation probability formalism). The latter 
implies dissolving of merging lines, which is the effect of non- 
ideal gas in a stellar atmosphere (see also Lanz & Hubeny, who 
applied theory of dissolving levels in hydrogen dominated DA 
atmospheres with carbon or iron in trace abundances). Such a 
nonideal gas effect has been included also in our tlusty195 
hydrogen atmosphere calculations, see Table 2 and Fig. 1 . 

Recent paper by Rohrmann et al. (2002) presents theory 
and derivation of the nonideal EOS for H and He gas, includ- 
ing their molecules. Rohrmann et al. (2002) concentrate on the 
theory of He core white dwarfs, which are much colder than 
WDs discussed in this paper. 

4. Pure iron and H-Fe stratified atmospheres 

We have computed and present here extensive set of pure iron 
model atmospheres and synthetic UBVRI broadband colors for 
models of 20000 < T eff < 100000 K and logg = 7, 8, 9, and 
10 (cgs). We have extended our computations to an extremely 
large logg. In fact we are not aware of the existence of any 
white dwarf star with log g > 10 but assume that iron cores of a 
given mass M would have much smaller radius R, and therefore 
log g, than common carbon-oxygen cores. Pure H model atmo- 
spheres discussed in the previous section were also extended to 
such a large surface gravity due to the same reason. 

Pure Fe model atmospheres were computed using nonideal 
EOS by Rogers (2002) and Los Alamos TOPS monochromatic 
opacities. Iron models were computed also on the grid of 155 
discrete standard optical depth points, starting from the ex- 
tremely high layer t\ = 10~ 10 and with 901 discrete frequen- 
cies. Computations were very difficult and sometimes unsta- 
ble due to complicated run of monochromatic opacities in Los 
Alamos tables and due to high gradients of opacities (with re- 
spect both to temperature and gas density) at higher densities of 
this exotic chemical composition. Moreover, pure iron model 
of T e /f = 100000 K and log g — 7 was not in hydrostatic equi- 
librium, since the model exceeded the Eddington limit in its 
deeper layers. In other words, the acceleration of gas, exerted 
by the integrated radiative flux and directed upwards, exceeded 
gravitational acceleration g in those layers. 

4.1. Physical conditions in Fe atmospheres 

Physical properties of pure Fe white dwarf atmospheres do not 
differ substantially from pure H atmospheres, at least for rela- 
tively high effective temperatures T e ff > 20000 K which were 
investigated in this paper. Table [3] presents concisely values 
of few most important thermodynamical variables across the 



sample pure Fe model atmosphere with T e ff = 20000 K and 
l°g 8 = 8.0, organized in the same way as it was described in 
Section 3.1. 

After each temperature iteration (and also after the final 
one), gas presures P g in the model are determined by integrat- 
ing the equation of hydrostatic equilibrium 

dPg_ _ g _ dPrad (4) 
dTstd (k v + cr v ) std dr std ' 

where all the monochromatic quantities were computed at the 
standard wavelength, A st d = 1500 A. One can note that gas 
pressure P g (T st d) is essentially determined by run of the stan- 
dard opacities, (k v + cr v ) st d- 

Column No. 6 in Table [3] shows that TOPS opacities (for 
1 gram of gas) are only « 4 times larger for iron than for hy- 
drogen. This implies that resulting gas pressures P g on a fixed 
r st d are few times smaller in pure Fe sample atmosphere than in 
pure H atmosphere of the same T e /f = 20000 K. Iron and hy- 
drogen atmospheres do not differ substantially from each other 
in the run of physical variables T, P g , and the density p. 

On the contrary to the pure H model, columns 4 and 5 of 
Table [3] show that the realistic nonideal Fe gas density proc 
(Rogers 2002) differs from the density pjo, computed by our 
atm21 code according to the ideal gas EOS. The difference of 
both densities reaches one order of magnitude in that model, 
which clearly demonstrates that the ideal EOS of iron would 
be of no use in WDs of T e ff = 20000 K. We have checked that 
iron gas of density p = 10~ 5 (cgs) behaves like an ideal gas in 
temperatures T as high as 10 6 K or higher. 

Column 7 in Table|3]presents geometrical depths of various 
standard optical depth levels r st d, assuming that the zero (refer- 
ence) level of geometrical height is fixed at the upper boundary 
of a model atmosphere. The depth of a fixed level r st d is given 

by 

Tad 

HGT = - f ; 7~ ; ■ ( 5 ) 

J p{K v + cr v ) st d 



The total thickness of a sample pure H model atmosphere of 
T e ff = 20000 K and \ogg = 8.0 equals to 7 kilometers. The 
corresponding pure Fe atmosphere is thinner, obviously due to 
a higher molecular weight of an iron gas. 

4.2. Spectra of Fe atmospheres 

Sample spectra of both pure iron and pure hydrogen model at- 
mospheres are plotted in Figs 2-7. They are arranged in order 
of increasing effective temperatures T e ff to demonstrate quanti- 
tatively evolution of both Fe and H spectra with increasing T e ff 
in the wide spectral region where the most of radiative energy 
is emitted. Effective temperatures T e ff chosen for this presen- 
tation range from 20000 K up to 100000 K and logg = 8.0 or 
9.0 in hotter models. 

Both series of model atmospheres were computed with the 
nonideal EOS from the OPAL Project to ensure consistency. 
One can see that pure Fe spectrum is entirely different from 
the H spectrum, both in visual (Fig. 3) and in far and extreme 
UV (Figs. 2 and 4-7). The Los Alamos TOPS monochromatic 
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Table 3. Model atmospheres of T e ff = 20000 K, and \ogg = 8.0 (cgs units), with pure hydrogen or pure iron composition. 
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opacities suitable to nonideal EOS allowed us to reproduce 
strongly broadened hydrogen Lyman and Balmer lines (Fig. 
3). Appearance of theoretical optical and near UV spectrum of 
the iron atmospheres is completely different and also exhibits 
strongly broadened iron features, but their identification is be- 
yond the scope of this paper. 

Fig. 3 shows that the sample pure Fe atmosphere of T e ff = 
20000 K is much brighter than pure H atmosphere in visual and 
also in far UV region beyond the Lyman edge (wavelengths 
A < 912 A). In near UV, for A > 912 A, pure Fe atmosphere is 
fainter that its pure H (i.e. the DA type) counterpart. Figs. 4-7 
demonstrate that Fe atmospheres are very bright for A < 912 
A also in higher T e ff. Only in the extreme UV, for A < 300 - 
400 A, Fe atmospheres get extremely faint apparently due to 
increasing iron monochromatic absorption. 

The latter effect, i.e. an increase of iron absorption in EUV 
and X-rays, is clearly seen in TOPS tables at any fixed tem- 
perature T and gas density p. In general, iron gas is much more 
opaque than hydrogen at any wavelength. This implies also that 
our pure Fe model atmospheres are geometrically thinner, as 
compared with pure H atmospheres of the same T e fj and log g, 
see also Table [5] 

4.3. Broadband UBVRI colors 

Theoretical broadband UBVRI colors of our pure iron models 
are displayed in Table 3. Comparison of Tables 1 and 3 shows 
that U - B and B - V colors of iron models are significantly 
lower than the same color indices of hydrogen model atmo- 
spheres. Therefore, a pure iron atmosphere of a given T e ff and 
log g appears as if it was much hotter than pure H atmosphere 
with the same parameters (an iron WD appears more blue to an 
observer). 

We have also computed model atmospheres of a more com- 
plex structure, in which pure iron gas is covered by a pure hy- 
drogen layer of an arbitrary thickness. Table 4 presents color 



Table 4. Broadband colors in stellar magnitudes for pure iron 
atmospheres. Models are computed with nonideal EOS and Los 
Alamos (TOPS) opacities. Units of T e fj and logg are the same 
as in Table 1 . 



z 


T eff 




U-B 


B-V 


V -R 


V-I 


26 


20 


7 


-1.330 


-0.284 


0.097 


-0.063 


26 


20 


8 


-1.188 


-0.191 


0.011 


-0.133 


26 


20 


9 


-1.110 


-0.111 


-0.047 


-0.185 


26 


20 
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26 
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26 


50 


9 
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-0.459 


26 


50 


10 


-1.440 


-0.235 


-0.188 


-0.418 


26 


70 


7 


-1.322 


-0.387 


-0.174 


-0.420 


26 


70 


8 


-1.311 


-0.396 


-0.188 


-0.443 


26 


70 


9 


-1.303 


-0.379 


-0.193 


-0.459 


26 


70 


10 


-1.313 


-0.340 


-0.197 


-0.467 


26 


100 


7 


-1.323 


-0.376 


-0.168 


-0.423 


26 


100 


8 


-1.315 


-0.358 


-0.133 


-0.321 


26 


100 


9 


-1.327 


-0.419 


-0.187 


-0.442 


26 


100 


10 


-1.321 


-0.425 


-0.195 


-0.461 



indices of such stratified atmospheres, in which we have arbi- 
trarily assumed that the transition depth between hydrogen and 
iron layers is located at the standard optical depth to = 10 3 . 
This assumption precisely determines the column mass of a hy- 
drogen layer in each model atmosphere 

TO 

col. mass = | (6) 

J (K v + cr v ) std 
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Table 5. Broadband colors in stellar magnitudes for stratified 
hydrogen/iron models with nonideal EOS and Los Alamos 
(TOPS) opacities. Units of T e ff and \ogg are the same as in 
Table 1. Transition depth is located at t = 10~ 3 . 
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-0.192 


H/Fe 


20 


10 


-1.098 
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Hydrogen column mass in our H-Fe stratified atmospheres 
weakly depends both on T e ff and logg of the models. For ex- 
ample, H column mass is equal to 9.6025 x 10~ 5 g cirT 2 in 
the coldest model of T e ff = 20000 K and \ogg = 8, and 
slowly rises to 3.0602 x 10~ 4 g cirT 2 for the hottest model of 
T e ff = 100000 K and log g = 9 (cf. sequence of sample models, 
Figs. 2-7). 

It would be very interesting to convert H column masses to 
the total mass of a hydrogen envelope covering the iron core of 
a white dwarf. Such a rescaling requires the knowledge of the 
radius R, whereas our models are parametrized by the gravity, 
g ~ M/R 2 . We have arbitrarily assumed the mass of a white 
dwarf to be M = M , and used M - R relations for iron core 
white dwarfs (Panei et al. 2000; Althaus 2003). Respectively, 
the total mass of a hydrogen envelope in two above sample 
models is lesser than 10 _19 M Q for T eff = 20000 K and logg = 
8, and is approximately equal to 2.6 x 10~ 19 M G for the hottest 
model of T eff = 100000 K and logg = 9. 

One can note that the estimated masses of a hydrogen 
layer are rather small, as compared with H masses arbitrar- 
ily assumed in published grids of theoretical M - R relations 
(M = 10~ 16 up to 10" 4 M o ). This is because H-Fe stratified 
models in the present paper require a H layer of the moder- 
ate optical depth t v < 1 in frequencies corresponding to the 
UBVRI filters, which do not cover completely the iron core. 

Figs. 8 and 10 show positions of various model atmo- 
spheres on the theoretical B - V vs. U - B diagram, whereas 
Figs. 9 and 1 1 display colors of the models on the B - V vs. 
V-I diagram. The Figures present colors of both pure and strat- 
ified atmospheres. For these models, we obtained U - B, B - V 



and V-I colors which depend on changes of T e ff and logg 
in a rather complex way (cf. Tables 3-4). One can clearly see 
that the colors of iron-containing models may overlap, causing 
them to be a non-unique function of the set of parameters T e ff, 
log g. Pure hydrogen models are also subject to this feature (at 
least on the B - V vs. V - I diagram), but the complexity for 
the iron models is much higher. This complicate dependence 
of the colors is due to a noisy structure of the TOPS grids of 
iron monochromatic opacities. Consequently, all these Figures 
present shaded areas, where atmospheres containing pure iron 
are located. We do not highlight lines of constant T e ff or log g 
for those models. One can see, that the region occupied by 
pure iron model atmospheres significantly differs from the re- 
gion corresponding to pure hydrogen models in all the analysed 
cases. 

The latter observation shows that it would be possible to 
distinguish pure iron WD atmospheres from pure hydrogen DA 
atmospheres with T e ff > 20000 K by means of the broadband 
Johnson UBV photometry, at least in principle. Our analysis 
and Figures include neither possible admixture of helium or 
other heavier elements in DA white dwarf atmospheres, nor 
the influence of their opacities on spectra and theoretical colors 
of mostly hydrogen white dwarf atmospheres. Admixture of 
heavier elements in hot WD atmospheres seems very common 
in real stars. For example, Barstow et al. (1993) showed that 
most of the hot DA white dwarfs with effective temperatures 
T e ff > 40000 K have some abundance of elements heavier than 
He in their atmospheres. This can explain their generally low 
X-ray luminosity, as it is found in the ROSAT all-sky catalog 
of white dwarfs. 

There exist many U - B and B — V color measurements 
available for isolated white dwarfs, which were published in 
the catalogue by McCook & Sion (1999). Therefore we have 
plotted colors for individual white dwarfs in both Figs. 8 and 
10 as asterisks, in order to search for possible white dwarfs with 
pure or stratified Fe atmospheres. Unfortunately, one can note 
that the observed stars generally do not fall into shaded areas of 
each Figure, but rather avoid them. Most of the observed U - B 
and B - V colors mostly overlap the area of pure hydrogen 
atmospheres. 

The above observation proves that the set of white dwarfs 
with color indices predicted for iron atmospheres is almost 
empty, except for a few cases of marginal significance. 

Figs. 9 and 1 1 do not include observed colors of real WDs. 
Measurements of near infrared luminosity / (and therefore V-I 
color) are extremely rare for hotter isolated white dwarfs. Also 
the newest catalogs of B - V and V-I of hot objects, which are 
known to us (Renzini 2001; Zoccali et al. 2001) do not distin- 
guish hydrogen DA and other types of white dwarfs. Existing 
surveys of V - I colors of white dwarfs concentrate on the 
search of infrared excess presumably in binary stars, and are 
not useful in this research. 
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Fig. 2. Pure iron and hydrogen LTE white dwarf synthetic spec- 
tra of T eff = 20000 K and logg = 8.0, displayed on wide 
range of wavelengths from near infrared, A = 10000 A, down 
to about 500 A. Solid line represents the spectrum of pure iron 
model and dotted line is the spectrum of pure H atmosphere. 
Both model spectra are computed using the EOS of nonideal 
gas and the Los Alamos TOPS monochromatic opacities. 



Fig. 4. Same as Fig. 2 for T eff = 30000 K and logg = 8.0. 
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Fig. 3. The LTE theoretical spectra of iron and hydrogen at- 
mospheres of Fig. 1 are displayed here on narrower range of 
wavelengths, with higher resolution, from near A = 3200 A up 
to 8000 A. Again, solid line represents the spectrum of a pure 
iron model and dotted line is the spectrum of a pure H atmo- 
sphere. 



5. Search for iron core white dwarf candidates 

5.1. Listing by Provencal et al. (1998) 

It is interesting to search for possible iron core white dwarf 
stars among the existing WD catalogs. The first such effort has 
been done by Provencal et al. (1998), who compared masses 
and radii of white dwarf stars with theoretical mass-radius re- 
lations for various WD core compositions. Their mass and ra- 




Fig.5. Same as Fig. 2 for T eff = 50000 K and logg = 9.0. 
Note, that the surface gravity is higher than in previous Figures. 




log wavelength (angstr) 



Fig. 6. Same as Fig. 2 for T eff = 70000 K and logg = 9.0. 
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Fig. 7. Same as Fig. 2 for T eff = 100000 K and logg = 9.0. 

dius determinations were based mostly on the Hipparcos par- 
allaxes, thus suggesting higher accuracy. Mentioned paper in- 
cludes also Procyon B as a candidate for iron core star, how- 
ever, this star has been classified as a normal white dwarf 
based on precise Hubble spectral measurements (Provencal et 
al. 2002). 

Based on the first quoted paper we have selected 7 white 
dwarfs, for which best determined masses and radii place 
them closest to the theoretical mass-radius relation for iron 
core white dwarf stars (cf. Fig. 3 in Provencal et al. 1998). 
Neglecting possible errors of M and R determinations we sug- 
gest that these stars are promising candidates for comparing 
them with our theoretical UBVRI color indices (see Table 5). 
Some of these stars were also indicated by Panei et al. (2000), 
who placed them on their theoretical mass-radius relation for 
iron core white dwarfs. 

Unfortunately, stars listed in our Table 6 are rather cold. 
The hottest of them, EG 50 (T eff = 21000 K) and GD 140 
{J e ff = 21700 K), are located practically at the border of the 
grid presented in this paper. Color indices and the position of 
both stars in our two-color diagrams places them far from the 
region occupied by pure iron atmospheres. Moreover, they are 
located in the region occupied by pure hydrogen atmospheres. 
Therefore we conclude that even if those stars have iron cores, 
they are covered by a thick hydrogen layer of a large optical 
depth. This conclusion is supported by the fact that McCook & 
Sion (1999) classified all the stars from Table 5 as pure hydro- 
gen DA white dwarfs. 
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Table 6. List of iron core candidate stars selected from the data 
by Provencal et al. (1998). Procyon B has been excluded al- 
ready from their data. 



EG 21 


0.554 


8.06 


0.0115 


-0.55 


0.05 


EG 50 


0.497 


8.10 


0.0104 


-0.91 


-0.08 


G 238-44 


0.417 


7.90 


0.0120 


-0.86 


-0.09 


G 181-B55 


0.50 


8.05 


0.0110 


-0.55 


0.17 


GD 140 


0.797 


8.48 


0.0085 


-0.95 


-0.06 


GD279 


0.411 


7.83 


0.0129 


-0.60 


0.10 


WD 2007-303 


0.433 


7.86 


0.0128 


-0.66 


0.07 



The above papers discuss mostly DA white dwarfs, i.e. 
those with hydrogen-rich atmospheres and therefore strong 
Balmer lines (cf. Figs. 2-3). Grids of theoretical UBVRI col- 
ors of Fe-abundant atmospheres given in Tables 3-4 are not 
useful here. Therefore we attempted to identify stars of excep- 
tionally high logg (and hence a small radius) in those papers. 
Unfortunately, several trial figures displaying M vs. logg re- 
lations have shown that there are no stars which exhibit large 
log g for a given mass as compared with the general trend. 

As an example, we plot two such relations in Fig. 12, which 
are based on the data from Finley et al. (1997), and Homeier et 
al. (1998). The former paper presents a table with mass and sur- 
face gravity determinations for 90 hot DA white dwarfs from 
the Extreme Ultraviolet Explorer (EUV E) all-sky survey. 




WD 



M/M Q \ogg R/Rq U-B B-V 



5.2. Other catalogs of WD masses and surface 
gravities 

In order to search for candidate iron core stars we have per- 
formed search of papers listing mass and surface gravity de- 
terminations of white dwarfs carried out in the recent surveys 
(Vennes et al. 1997; Finley et al. 1997; Homeier et al. 1998; 
Napiwotzki, Green & Saffer 1999; Vennes 1999; Bergeron et 
al. 2001). 
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Johnson UBV color grids for model atmospheres: 

pure iron (filled / shaded area); 
pure hydrogen (right-upper side). 

Stars represent McCook & Sion 
(1999) catalog data for white dwarfs. 
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Fig. 8. Plot of U - B vs. B - V for 
hot white dwarf stars. Shaded area de- 
notes region occupied by colors of pure 
iron model atmospheres computed in 
this paper. Numbers put in this area 
denote effective temperatures T e ff (in 
thousands K) and logg, respectively. 
Colors of pure hydrogen models (non- 
ideal EOS) are represented by the trans- 
parent grid of sections. Asterisks de- 
note colors of real white dwarfs, as 
taken from McCook & Sion (1999). 
One can see that only few real white 
dwarfs fall into the shaded area. 
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Fig. 9. Plot of B - V vs. V - 1 for hot 
white dwarf stars. Also here, shaded 
area and the adjacent numbers corre- 
spond to the colors of pure Fe model 
atmospheres. Transparent grid of sec- 
tions correspond to pure H models 
(nonideal EOS). Unfortunately, there 
are no observed color indices available. 
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Johnson UBV color grids for model atmospheres: 

stratified H / Fe (filled / shaded area 
pure hydrogen (right-upper side). 

Stars represent McCook & Sion 
(1999) catalog data for white dwarfs 
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Fig. 10. Plot of U - B vs. B - V for 

hot stratified Fe/H model atmospheres. 
In these models, pure Fe atmosphere 
is covered by a pure H layer, reach- 
ing from the standard optical depth 
t = 1(T 10 down to 10~ 3 in all cases. 
Designations of T e ff and log g have the 
same meaning as in previous Figures. 
Pure H models are also represented in 
the same way. 
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Fig. 12. Catalog of M/M vs. log g determinations for two sets 
of DA white dwarfs. Dots denote stars observed by EUVE 
satellite by Finley, Koester, & Basri (1997). Crosses represent 
data from the Hamburg Quasar Survey (Homeier et al. 1998). 
All the investigated stars are located in well defined strip cor- 
responding to C-O white dwarf cores. No single star with too 
high log g is seen in the Figure. 



spheres. We showed that effects of nonideal EOS for pure iron 
model atmospheres are substantial in the investigated models, 
contrary to pure H WD atmospheres. 

4. We have also computed UBVRI colors of the models 
and determined the area on the B - V vs. U - B and B - V 
vs. V - I planes, occupied by both pure Fe and pure H model 
atmospheres of WD stars. 

The importance of our paper is that we have determined 
the area on two-color diagrams in which colors of the hypo- 
thetical pure iron or stratified H/Fe atmospheres should be lo- 
cated. Colors of the above hypothetical atmospheres are dis- 
tinctly different from colors of pure hydrogen DA WSs with 
T e ff > 20000 K. We are aware, however, that our paper does 
not predict colors of iron core WDs in case, if the core is cov- 
ered by a hydrogen layer of significant optical thickness. 
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Surface gravity \ogg and T e ff are usually determined by 
fitting of the observed Balmer line profiles and therefore this 
method makes direct to observable quantities. One should nec- 
essarily note, however, that in all the above papers masses of 
DA white dwarfs were obtained using the theoretical mass- 
radius relations for carbon-core stars with hydrogen or helium 
envelope (see Wood 1995). We are aware that if masses were 
determined with the most new mass-radius relation for iron- 
core WDs by Panei et al. (2000), then the above Figure would 
change significantly. 

6. Summary 

The most important areas of concern of our paper can be sum- 
marized as follows. 

1. We have computed and presented here LTE model at- 
mosphere computations of hot white dwarfs, which atmo- 
spheres consist either of pure iron, stratified iron/hydrogen lay- 
ers, or pure hydrogen. Model atmospheres are in both radiative 
and hydrostatic equilibrium. Computations use exclusively the 
equation of state for nonideal Fe and H gases, obtained by the 
OPAL project, and the extensive tables of the TOPS monochro- 
matic opacities from the Los Alamos National Laboratory. The 
latter opacities are fully compatible with the nonideal EOS for 
both gases. 

2. This paper presents calibration of colors U - B, B - V, 
V - R and V - I of the broadband Johnson UBVRI photome- 
try, for pure Fe, pure H and stratified Fe/H atmospheres. These 
colors were computed from theoretical spectra of our model 
atmospheres with a nonideal EOS. 

3. The paper investigates the influence of nonideal equa- 
tion of state on the photometric UBVRI colors of the broad- 
band Johnson system, both in pure H and pure Fe model atmo- 
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